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It is well-known that effective field theories with four-fermion interaction (the so-called Nambu 
- Jona-Lasinio (NJL) models), which incorporate the phenomenon of dynamical chiral symmetry 
^ \ breaking, are quite useful in describing low-energy hadronic processes (see e.g. [H HI and references 
therein). Since the NJL model displays the same symmetries as QCD, it can be successfully used 
for simulating some of the QCD ground state properties under the influence of external conditions 
such as temperature, baryonic chemical potential, or even curved space-time etc. |E IE IE IE El I n 
particular, the role of the NJL approach significantly increases, when numerical lattice calculations 
^ ■ are not admissible in QCD, like at nonzero baryon density and in the presence of external gauge 
^p" 1 - fields |EIH- I n this way, it was observed in the framework of a (2+1) -dimensional NJL model that 
Ph! an arbitrary small external magnetic field induces the spontaneous chiral symmetry breaking (CSB) 
even at arbitrary weak interaction between fermions Later, it was shown that this phenomenon 
(the so called magnetic catalysis effect) dis pla ys a universal character and can be explained basing 
on the dimensional reduction mechanism 0, ^[ (for a modern status of the magnetic catalysis 
effect and its applications in different branches of physics, see the reviews (l2|). 

As an effective theory for low energy QCD, the NJL model does not contain any dynamical gluon 
fields. However, such a nonperturbative feature of the real QCD vacuum, as the nonzero gluon 
condensate (F^ l/ F atJ-u ) can be mimicked by external chromomagnetic fields. In particular, for a QCD- 
motivated NJL model with gluon condensate (i.e. in the presence of an external chromomagnetic 
field) and finite temperature, it was shown that a weak gluon condensate plays a stabilizing role 
for the behavior of the constituent quark mass, the quark condensate, meson masses and coupling 
constants for varying temperature 13]. Then, in a series of papers, devoted to the NJL model with 
the gluon condensate, it was shown that an external chromomagnetic field, similar to the ordinary 
magnetic field, serves as a catalyzing factor in the fermion mass generation and dynamical breaking 
of chiral symmetry as well 14]. The basis for this phenomenon is also the effective reduction of the 
space dimensionality in the presence of strong external chromomagnetic fields, and this does not 
depend on the particular form of the field configurations [15[ . 

At present time, it is well-established (see, e.g., []j|]) that the gluon condensate is a very slowly 
decreasing function of the baryon density (baryon chemical potential). So in cold quark matter, it is 
a nonzero quantity even at sufficiently large baryon densities, which are expected to exist inside the 
neutron star cores. Evidently, the consideration of the gluon condensate may change significantly 
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the generally accepted picture of physical processes in the baryon matter. For instance, according to 
the modern point of view, the color superconducting (CS C) p henomenon can be realized in neutron 
star cores (see the reviews pj). However, according to 17], the critical parameters of this phase 
transition strongly depend on the value of the gluon condensate (in this case, only its chromomagnetic 
component survives). Moreover, it turns out that even for a rather weak quark coupling, different 
external chromomagnetic field configurations induce the CSC phenomenon ^j]. Quite recently it 
was found that some combinations of external chromomagnetic and ordinary magnetic fields can 
penetrate into a bulk of the CSC medium and modify its ground state, producing a new type of 
color superconductivity [l^]. Finally, note that in the dense quark matter, a superstrong magnetic 
field is originated due to the presence of the gluon condensate j^j. In conclusion, we see that there 
exist several new physical effects that are intrinsically connected with the gluon condensate (external 
chromomagnetic fields). 

In the present paper, we study the influence of an external chromomagnetic field on the pion 
condensation phenomenon. This phase transition can also occur in the dense baryon matter, although 
isotopic asymmetry, with different densities of up and down quarks, is needed for this process to 
take place. This type of quark matter was already investi gate d in the framework of NJL type of 
models, both with and without pion condensation jU HH The main purpose of the 

present paper is to investigate, in the framework of a NJL model, the behavior of the quark and pion 
condensation in the presence of an external chromomagnetic field, modelling the gluon condensate. 
In particular, we will show that, even for a weak coupling of quarks, the pion condensation effect is 
induced by an external chromomagnetic field and is related to an effective dimensional reduction. 
The latter effect leads to a nonanalytic logarithmic dependence of the quark and pion condensates 
on the field strength in the strong field limit. 

II. QUARK AND PION CONDENSATES IN EXTERNAL FIELDS 

For a system, composed of two flavored quarks, there exists the relation uq = nj 3 + between 
the electric charge density uq, the baryon charge density ub and the density nj 3 = \qr^q of the 
third isospin component 1% = T3/2. Since these quantities are linearly dependent, we will study in 
the following the isotopic asymmetry (which means that different species of quarks have different 
densities) in dense quark matter with a nonzero baryon chemical potential fis- in this case, the 
chemical potentials fiB and fij (where fij is the isospin chemical potential) are independent parame- 
ters. Note that in another possible case, i.e., \ib 7^ 0, \iq 7^ (the last quantity is the electric charge 
chemical potential), they are no more independent, but related through the electric charge neutrality 
condition nq = 0, if one assumes that quarks are in a weak-equilibrium with electrons, and the whole 
system is electrically neutral. We restrict ourselves here to the first possibility, i.e. suppose that hb 
and hi are independent quantities. In our investigations we also suppose that \ib is a rather small 
quantity. In this well as at sufficiently high flavor asymmetry, the color superconductivity 

effects may be neglected in favour of the normal quark matter or pion superfluidity effects (see the 
recent discussions in 0, 0|)- 

A. General definitions 

Let us consider a NJL model of flavored and colored quarks q ia {i = 1, . . . , Nf, a = 1, . . . , N c ) 
with Nf — 2, N c — 3 as numbers of flavors and colors, respectively (for convenience, corresponding 
indices are sometimes suppressed in what follows), moving in an external chromomagnetic field. 
The underlying quark Lagrangian is chosen to contain four-quark interaction terms responsible 
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for spontaneous breaking of chiral and flavor symmetries. Hence, two types of condensates might 
characterize the ground state of the model: the quark condensate (qq) (spontaneous breaking of 
chiral symmetry), and the pion condensate (qj5T±q) (spontaneous breaking of parity and isotopic 
symmetry). In particular, we consider a Lagrangian which describes dense quark matter with an 
isotopical asymmetry, and we neglect diquark interactions and hence the possible formation of a 
diquark condensate. Upon performing the usual bosonization procedure j2f|, [l[ and introducing 
meson fields a, n, the four-quark terms are replaced by Yukawa interactions of quarks with these 
fields, and the Lagrangian takes the following form (our notations refer to four-dimensional Euclidean 
space with it = £4 r ): 

£ = -q(i^ u S7 u 4- inj Q 4- a + m 4- i^Tii + ifi'r 3 j°)q - -^(cr 2 + 7? 2 ). (1) 

Here \i = (fi u + /id)/2 is the chemical potential averaged over flavors 2 , // = (fi u — fid)/ '2 is their 
difference, and G is the (positive) four-quark coupling constant. Furthermore, V M = <9 M — igA a ^\ a /2 
is the covariant derivative of quark fields in the background field = d^A^ — d u A a ^ + gfab c A b „A° 
determined by the potentials A a (a = 1, 8), and A a /2 are the generators of the color S77 c (3) group. 
Finally, r = (r 1 , r 2 , r 3 ) are Pauli matrices in the flavor space. 

Evidently, the Lagrangian {!]) with nonvanishing current quark masses, m 7^ 0, is invariant under 
the baryon £7b(1) symmetry and the parity transformation P. Moreover, without the //-term, it is 
also invariant under the isospin SUi{2) group which is reduced to Ui 3 (l) at fjf 7^ 0. At m = and 
fi' 7^ the symmetry group of the initial model is Ub(1) x Ui 3 l(X) x [/^(l) x P (here the subscripts 
L,R mean that the corresponding group acts only on the left, right handed spinors, respectively). 
It is very convenient to present this symmetry as £/b(1) x Uj 3 (1) x Uai 3 {^) x P, where Uj 3 (l) is 
the isospin subgroup and Uai 3 {]) is the axial isospin subgroup. Quarks are tranformed under these 
groups in the following ways q — ► exp(iar 3 )g and g — > exp(iccy 5 r 3 )q, respectively. In the case of 
m = the phase structure of the model ((IJ is defined by the competition of only two condensates. 
One of them is the quark condensate (qq), and the other is the pion condensate (qjsTxq). If the 
ground state of the model is characterized by (qq) 7^ and (q^T\q) = 0, then the axial symmetry 
Uai 3 (1) of the model is spontaneously broken down, but isospin symmetry Ui 3 (l) and parity P remain 
intact. However, if (q^Tiq) 7^ and (qq) = 0, then only parity P and the isospin symmetry are 
spontaneously broken down and the pion condensed phase is realized in the model. As a consequence, 
in the last case we have (qr 3 q) 7^ 0, i.e., a nonzero difference in densities of up and down quarks arises 
(isotopic asymmetry of the ground state). Generally speaking, the inverse is not true, i.e., if there 
is an isotopic asymmetry of quark matter, this does not necessarily mean that pion condensation 
phenomenon does occur. 

In order to investigate the possible generation of quark and pion condensates in the framework of 
the initial model (JT|). let us introduce the partition function Z of the system 



exp We = / dqdqdadiii exp 



d xC 



(2) 



1 We consider 7— matrices in the 4— dimensional Euclidean space with the metric tensor = diag(— 1, — 1, — 1, — 1), and the relation 
between the Euclidean and Minkowski time = * x (m) : = ^(E) = 7(V)' ^ n wna * follows we denote the Euclidean 
Dirac matrices as 7 M , suppressing the subscript (E). They have the following basic properties 7^" = — 7 M , {7^,7^} = — 2<5 M „. The 
charge conjugation operation for Dirac spinors is defined as Vc(^) = C (ip{ x ) + ) t with C^C~ 1 = — 7^. We choose the standard 
representation for the Dirac matrices (see 1271 ). The 75 has the following properties: {7 M ,7s} = 0, 7^" = 7g = 75- Hence, one finds 
for the charge-conjugation matrix: C = 7°7 2 , C + = C~ 1 = C* = — C. 

2 It is equal to one third of the baryon chemical potential: fi = (j.b/3. Since the generator ^3 of the third component of isospin is equal 
to T3/2, the quantity fi' in is half the isospin chemical potential, fi' = /z//2. 



with We being the Euclidean effective action. Here, the meson fields can be decomposed as follows 

a = ct (0) +5a, 7f = tT (0) + 5n, 

where and 7f(°) are the condensate fields, which are determined by the minimum of the effective 
action, 

Notice that the above expectation values are expressed through functional integrals over the quark 
and meson fields, where the latter are calculated by using the saddle point approximation (in what 
follows, we shall denote these mean values again by a and tt, respectively). Thus, instead of Eq. (J2J), 
we shall deal with the following functional integral over fermion fields: 



exp We = / dqdq exp 



d xC 



(3) 



where now a and tt in the Lagrangian (0) are understood as constant condensates. In the mean field 
approximation, where field fluctuations 5a and 6n are neglected, they are given by the following gap 
equations 

Assume that the only nonvanishing components of the background gauge field potential are for 
a = 1,2,3, while others are equal to zero, 3 

^ 0, when a = 1, 2, 3, and A* = 0, when a = 4, . . . , 8. 

This implies that only quarks of two colors a = 1, 2 do interact with the background field A a . As a 
result, the integration over quark degrees of freedom in the partition function © is greatly simplified, 
and we have 

Z = Det(i)(z7(9 + M) ■ Det( 2 )(ryV + M), (5) 

where M. = m + a + 17 7rr + i/^7o + ^'^37°, and indices (1) and (2) mean that determinants 
are calculated in the one-dimensional (with color a = 3) and in the two-dimensional (with colors 
a = 1, 2) subspaces of the color group, respectively. Assume now that the background field is 
constant and homogeneous, = const. Then the Dirac equations 

{i^d + M)ip = 0, {i^V + M)ip = 



Of course, if the external constant homogeneous gauge field fills the whole space, then the color as well as the rotational symmetries 
of the system are broken. This might be considered as an artifact, since in reality there exist both color and rotational invariance (we 
ignore the color superconductivity effects). In our case, in order to deal with a physically acceptable ground state, i.e. without color 
symmetry breaking, one can interpret it as a space split into an infinite number of domains with macroscopic extension. Inside each 
such domain there exist a homogeneous background chromomagnetic field, which generates a nonzero gluon condensate {F^ v F a ^ v ) . 
Moreover, the direction of the gauge field varies from domain to domain, so the averaging over all domains results in a zero background 
chromomagnetic field. Therefore, color as well as rotational symmetries are not broken. Strictly speaking, our following calculations 
refer to some given macroscopic domain. The obtained results turn out to depend on color and rotational invariant quantities only, and 
are independent of the particular domain. 



have stationary solutions with the energy spectrum e for quarks of flavor i and color a = 1,2,3 with 
quantum numbers k moving in the constant background field F® u (a = 1, 2, 3). In this case, we arrive 
at the following Euclidean effective action: 
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Here, r is the imaginary time interval, the sum is over the signs A = ±1 of the chemical potential 
//, the signs k — ±1 of the chemical potential corresponding to charge conjugate contributions 
of quarks, color indices a = 1,2, 3, and also over quantum numbers k of quarks, with a = 3 for free 
quarks and the spectrum 



e = Eftx = y (y / (fr + m) 2 +p 2 + vr| + A// 
and with a = 1, 2 for quarks with the spectrum 



+ 7r 2 + 7r 2 , (A = ±i; 



(7) 



; = e k;X = y (y / (a + m) 2 + n2 + 7r| + A//)' + vr 2 + vr 2 , (A = ±1), (8) 

moving in the background color field F^ u (a = 1,2,3). In the above formula, II 2 . stands for the 

eigenvalues of the squared Dirac operator —(7V) 2 with V = d — igA a \ a /2. 

In the case of finite temperature T = 1/(3 > 0, the thermodynamic potential Q = —We/ (PL 3 ) is 
obtained after substituting p 4 — > ^-(l + |), I = 0, ±1, ±2, 



/3 

i=+oo 



75EE E 10 



A, ft fc,a Z=— 00 



2vr(/ + 1/2) 





+ (e — «/u) 2 



+ 



a^vr 2 



4G 



C, 



(9) 



where we introduced a subtraction constant C in such a way that at a = 7? = we have f2 = 0. 
Next, let us consider the proper time representation 
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( t 2 + tT 2 
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(10) 



1/AJ 



where A s is an ultraviolet cutoff (A s ^> a, \tt\). The temperature dependent contribution can be 
further transformed with the help of the formula 



^2 exp[-s(27r///3 + x) 



[l + 2^exp(-^-)cos(x/3/)], 



2 A /7TS 



l = — 00 i=l 

where in our case x = j|. Then, calculating the quark condensate 
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and combining the result of ()12j) and (jlj), we obtain the following gap equation: 



(12) 
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Here, the first term in the square brackets describes the T = contribution (it corresponds to 
the result of integration over p, in the initial equation ©), while the second term is the finite 
temperature contribution, T / 0. The 7f-condensate can be obtained in a similar way. For tt 3 we 
have 
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G 
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£ 3 A 



X,k k,a 



1/A§ 



1 + 2JJ, 
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7Ti* 



7TS 
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■ (14) 



By comparing this expression with that for a, we conclude that for nonvanishing quark mass m/fl 
this condensate vanishes, = 0. For tti (putting 7T2 = by consideration of the symmetry of the 
problem) we have 
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l + 2^(-l) i e 4 S 

Z=l 



(e — K/i) _ 



s (e 



(15) 



It follows from ()13j) and ()15j) that at m = and //' 7^ there exist only two different solutions of 
this system of gap equations (except for the trivial one with a — 7i"i = 0), i.e., a) a = 0, n 1 7^ and 
b) a 7^ 0, 71"! = 0. Thus we have to find out which of these solutions provide the global minimum of 
the thermodynamic potential © with fjf 7^ 0. 

In the limit of a vanishing external field (F° u = 0), with = 0, 7C2 = 0, we have for the particle 
spectrum 



and for the sum over quantum states 



(16) 



(2^ 



for 3 color states (a = 1,2,3), 2 spin states, and 2 values of k = ±1. Considering now the case of 
vanishing temperature T = 0, we shall omit the cutoff in the lower limit, 1/A 2 — > 0, and replace it 
by the corresponding cutoff in the momentum integration p 1 < A 2 . It can be easily seen that the 
two cutoff parameters A p and A s are related as A 2 = 2A 2 . Indeed, integration in (0) gives: 

~ J dp, J d s plog(e 2 + pi) = -J^At + 0(A 2 ^ 



57T 



p 2 <A 2 



with the momentum cutoff p 1 < A 2 and integration in (fTUj) gives 



00 



(27T) 



d 4 p 



dz 



exp(—ize 2 — izp 2 ) 



32tt 



-M + 0{A] 



1/A 2 



with the proper time cutoff z > In what follows, we shall denote the A 2 cutoff by A 2 , and we 
shall use it throughout, substituting A 2 = 2A 2 . Now, we can perform the proper-time integration in 
(fTSJl with the help of the formula 



K=±l 



ds (e - Kfi) c _ s ( £ _ KU )2 
y/s e 



\£ — fji\ 



7T- 



{e-fl) 



(17) 
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As a result, one obtains for the quark and pion condensates in the vanishing color field background 
the formulas coinciding with the result of 25], 

2G = 6 < ff + '">E J ^f— jWTWTf ' (18) 

711 _ C_ / A g ( g - /f) / 19 N 



2G ^7 (2vr) 3 e ' 

In what follows, we shall analyze the special case of a constant Abelian chromomagnetic field 

A" = SSS^H. (20) 

The spectrum IT| of the Dirac operator —(7V) 2 has then six branches, two of which correspond to 
quarks that do not interact with the chromomagnetic field (a = 3) 

n 2 , 2 = /, (21) 

and the other four correspond to two color degrees of freedom of quarks with "charges" ±g/2 
interacting with the external field. The Il| spectrum of quarks is now given by (a = 1, 2) 

^,4,5,6 = gH{n+ l - + ^)+ P l (22) 

where ( = ±1 is the spin projection on the external field direction, p% is the longitudinal component 
of the quark momentum (—00 < p 3 < 00), 

pl = gH (n+~) (23) 

is the transversal component squared of the quark momentum, and n = 0,1,2,... is the Landau 
quantum number. 

The form of the spectrum is essential for the quark and pion condensate formation. Using the 
above expressions for energy spectra, we shall next study the quark and pion condensates in the 
strong field limit. 



B. Asymptotic estimates for strong fields 



In this section, we consider the special cases of the above configurations of background fields 
in the strong field limit. Our goal is here to demonstrate that the field is a catalyzing agent for 
dynamical symmetry breaking, leading to a possible creation of corresponding condensates. The 
external fields are assumed to be strong as compared to the values of quark (qq) and pion (q^Tiq) 
condensates that may be rather small. In this sense, even the expected values of fields simulating 
the presence of a gluon condensate, which we take to be of the order of gH = 0.4 — 0.6 GeV 2 , may be 
considered as strong (these values of the external chromomagnetic field gH correspond to the QCD 
gluon condensate at zero temperature and zero baryon density [28]). As for the other parameters, 
we may take their values as in 0,H|, i.e., A = 0.65 GeV, G = 5.01 GeV -2 . 

Consider now the special choice of parameters: fi = 7r 2 = 7r 3 = m = 0. In this simple massless 
case, though with 7Ti 7^ 0, a 7^ 0, // 7^ 0, the thermodynamic potential (fTUj) takes the form 
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Next, consider the zero temperature case, T = 0. Then the second term in the square bracket in 
the above expression vanishes and we are left with the first term equal to unity. The above equation 
describes the effective potential V e R = fl\T=o 



V eS (a, V 
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2vW S n (2tt) 3 
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The first term in the square bracket stands for quarks with a = 3, when we have 



AG 



C. (25) 



4,A = (V / ^+^+A/i') 2 + vr 1 2 . 



(26) 



and 4 stands for two values of k — ±1 and two values of spin projection ( = ±1. The second term 
stands for quarks with a = 1,2, where we have 



fc fc,A 



(^x 2 + n 2 + A//) 2 



7Ti 



(27) 



Let us first consider the symmetric case, when the chemical potentials for the w-quark and the 
rf-quark are equal, // = (fi u — Hd)j1 = 0. Then, the energy of the free quark with a = 3 becomes 



— Q i 2 , 2 



-a , 2 



where a\ = a 2 + 

Consider now the case, when the Abelian-like background field is strong enough, g H = 0(A 2 ), 
but g H < 2A 2 . The momentum squared fl 2 , is given by ([22)1 . and summation over quark quantum 
numbers in the chromomagnetic field gives 



L 3 

X k,ct 



2 fEI> 

A ra=0 



5: 



nO 



dp 3 
2ir ' 



(28) 



Now, the momentum integral in the first term in (|23j) gives 



d 3 p -sel x 
e P' A 



(27T) 



7T e- S(7 i 
8tt 2 s 3 / 2 



In the second term in (|25j) for quarks with a = 1,2, due to the strong field condition, we take 
only the contribution of the term with n = in the sum over Landau quantum numbers, while the 
contribution of large quantum numbers is described similar to the contribution of the free quark 
term, and thus we must take the free term three times. As a result, after integration in the second 
term over the p3-component of the quark momentum we have 



4vr 2 



, -srr 2 3 
ds e s<T i — + 



sa 2 9H 
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The cutoff momentum A as well as the background field gH are large A ^> ai, gH ^> a\ and, hence 
we can approximately integrate over s in the first and the second terms according to the following 
formulas: 



dx 



-X 



1r 2 
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- ^log^n 



+ const, where x Tl 



2A 2 
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(29) 
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dx 



-x 



h logx min + const, where x min = — jz <C 1. 



(30) 



As a result we find the effective potential 
1 



47T 2 



3^ 



1 /2A 2 \ 2 2A 2 1, 2A 2 ^ 
+ - log — + d 



or 
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f + log ^ + C2 )} + il +c - 



< 31 > 

where Ci, C2 are certain numerical constants. Now, let us find the minimum do of the thermody- 
namic potential Ijfllj) . Then, the minimum o"o is described by the solution of the equation, where we 
neglect the terms that do not contain large parameters A 2 or gH log (giJ/cr 2 ) , 



cff 



-6A 2 -3a 2 log^ 



2# 



2 



<9<x 2 |CTl=ff0 4tt 2 
Within the approximation of strong background field 

gH\og{gH/al) » a\ log(2A 2 /a 2 ) 
the above equation simplifies to the following form: 



6GA 2 



71 z 



G— rloe^r 



(7n 



with the solution 



2tt 



exp 



2vr 2 (l-^) 
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where the effective coupling constant, 



6GA 2 



7T 



(32) 



(33) 



(34) 



(35) 



can be arbitrary small g < 1, contrary to the zero magnetic field case H = 0, when it should be g > 1 
in order that the nontrivial solution for o"o exists. This is a very interesting result. Indeed, we see 
that at y! = and arbitrary small attraction between quarks (g < 1) the external chromomagnetic 
field catalyzes the dynamical chiral symmetry breaking and quarks acquire a nonzero mass, which 
is equal to <j (|34|). In contrast, if H = and # < 1, the quarks remain massless and symmetry is 
unbroken. After removing the cutoff with the use of (|3*3*|) . the effective potential (|31j) takes the form 



gHa\ 



M 

V ^0 



(36) 



which is just the effective potential of the two-dimensional Gross- Neveu model (see, e.g., |29|)- This 
result shows that the chromomagnetic catalysis effect is accompanied by the effective reduction of 
the space-time dimensionality. 

Next, consider the case of flavor asymmetric quark matter with y! ^ 0. For simplicity, /i is taken 
to be zero as before. In this case (recall that m — 0), there are two nontrivial solutions of the 
gap equations (fT3*|) and (|T3j) (see the remark after (jl5|0 . which are the points on the a- or 7Ti-axes. 
So, in order to find the global minimum point, it is sufficient (and very convenient) to reduce the 
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investigation of the effective potential V e s (a, tti) (|2l)|) as a function of two variables o and 7i"i to two 
particular cases. First, we shall study V e ^ as a function of n 1 only with a = 0, 7i"i 7^ 0, then as a 
function of o with <r 7^ 0, 7i"i = 0. When comparing these two particular cases, one can obtain the 
genuine global minimum point of the effective potential V^(a,7Ti). To simplify our calculation, we 
assume that /1' and gH are large, /1' 2 = 0(A 2 ) and g H = 0(A 2 ). 

a) The case a = 0, n\ 7^ 0. 

In this case the first term in (|25j) . 



5'2 



(2tt 



(37) 



can be expressed through the following approximation of the integral (with s/i' 2 ^> 1): 



E 



^.e- s 4,A 



v^e- S7r i 



(2vr) 3 



4tt 2 



2 ^' 2 + T J -T- 
2s J Vs 



1 



1 



(3? 



Integration over s should be approximately performed separately for the first and second terms in 
(138)1 with different lower limits 



4vr 2 



4 



dsfi 



ts — snf 



+3 



ds 



-S7T7 



2A2" 



/ 



Here, we also added the contribution of quarks with a = 1, 2 with high Landau quantum numbers 
n > 1. 

In the second term in (|25|). 



1 



2^ 



(is 1 x ^ ^ 

~J3 1^ 



S 2 



e s£fc >\ 



l/A? 



k fe,A,o=l,2 



with gi/ = 0(A 2 ) and /x' = O(A), but gH < /i' 2 < 2A 2 , we should estimate the contribution of the 
two main terms, n = and n — 1, in the sum over n. The n = contribution is expressed through 
the integral 



<#>3 
(27T) 



exp[-s(7r 2 + (|p 3 | + A^) )] 



i e - S 7r 2 

7T 



(39) 



and summation over the remaining quantum numbers is made according to (j28j) with n — 0. For 
the n—1 term, we have a new expression with integration of the following exponential term: 



E 



^expHK 2 + (yW^ + V) 2 ) 



(40) 



In order to estimate the n = 1 contribution, we write P3 = ^3 + ip, where ip is a small deviation of 
£>3 from the minimum of the exponential quantity \J p\ + <?H = // for A = — 1. Then the integral in 
(J4T?|) for n = 1 is estimated as 
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Recall that this estimation is justified under the above assumption, i.e., gH = 0(A 2 ) and fi' = 0(A), 
but gH < fi' 2 < 2A 2 . Finally, summing up the n = and n = 1 contributions ()39j) and (}4~Tj) . we find 



K ( ff 2 Vi) 




In the second term in the above bracket, we took into account that the state with n — 1 is two- fold 
degenerate in the spin variable. Summing up all the contributions, and taking into account the 
integrals (}2T?j) and (jHU|) . we have 



Kff(7Tl) 



4tt 2 1 2 



2A 2 

7Ti 



2A 2 



7Ti 



log 



*1 

2A 2 



+ 



,4 



/'' 



m' 2 - gH 



+ 



I/'' (log^-0 2 



gHrr 2 



1 + 



/'' 



log | + C3 )} + ^ + a (42) 



The minimum point 7r of this function obeys the following stationarity equation dV e fi / dn 2 ^^^ = 0, 



which gives 



6GA 2 4/i ,2 G 



log^-^Gll 



g H . 



IT" 



7T 



/2 



2tt 2 



Vi^H) l0g 9H' 



(43) 



where the approximation of strong background field, gH\og(gH/nf) 7r 2 log(2A 2 /7r 2 ), has been 
used, similar to (J3~2~j) . After substituting the above equation into (H2J), the effective potential takes 
the form 



n 

Air 2 



4^ 2 + ^(l+ - " 

2 \ ^^gH 



log -1-1 



(44) 



which again, like (|36|). resembles the effective potential in the two-dimensional Gross - Neveu model 
. Note also that the minimum tcq of the function ()42|1. under our assumption of a strong color 
background field, exists even for weak coupling of quarks, i.e., at g < 1, while with zero background 
field, it can exist only if g > 1. Next, let us study the behaviour of the effective potential (}25j) as a 
function of a, when 7Ti = 0. 
b) The case a ^ 0, ix 1 = 0. 

Here again we use the above assumption gH = 0(A 2 ) and ^ = O(A), but gH < ^i' 2 < 2A 2 . The 
term in the effective potential with a = 3 is determined by the expression 



(i) 



off 



ds 

J si 

VA2 



4 £ 



(2vr 



(45) 



where = \J a 2 + p 2 + A//. The main contribution to the integral over p, for large //, is given by 
large p near p = a/ /i' 2 — c 2 . Then the lower limit in the integral over s should be taken as , 
and we obtain, with p = p + ip, 



V. 



(i> 



P 2 



cff 



7T 



5/2 



00 00 



(46) 
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Then, the first term, where the contribution of quarks with a = 1,2,3 in the region of very large 
1 

m' 2 ' gH 



values of s ^> , 4n is also added, is estimated as 



Wff) = -?(^-^)V + 



7T 



8tt 2 



2A 2 \ 2 A 2 , a 2 ' 



(47) 



Consider now the contribution of the term with a = 1, 2 in the region of comparatively small values 
of s close to the threshold 1/gH, 



V r 



cS2 



CO 

4^ / f E^E(2"M / ^expMy^T^T^n + V) 2 ]. (48) 



i/gH 



Since we assume that \s! 2 and (7 if are of the same order of magnitude and large, // 2 = 0(A 2 ) gH = 
0(A 2 ), two terms, n = and n — 1, will make the main contribution to the sum over n. The n = 
term is determined by the integral, approximately equal to 



A* 



/ ^exph^^/^i-Zi') 2 ]-- _ ^ r 
7 27r v 2 v /7ri A /Ai /2 - c 2 

The n = 1 term gives 

2 E / ^exph^^+^ + ^ + A/i') 2 ] 



where we had introduced a factor two due to the two-fold degeneracy of the n — 1 Landau level in 
the spin variable. The above two terms with n = and n — 1 should be summed up to give 



KflQ (cr) = 



8tt 2 



00 oc 

1 r ds 2 > 

72 ^ 2 7 s 2 + ^fJ^gH - a 2 J 



\ 



ds 



ff#A*' 
8tt 2 



// 



/2 



- a 2 + 2 V 7 / 2 -gH-a 2 ^ . 



(49) 



The effective potential after summing the contributions V e s\ and V e w, is finally estimated as follows: 



V eS (a) 



//' 



7T 



(fi' 2 -a 2 )i + 



3^ 
8tt 2 



2A 2 



-V- l0g 2A^ + Cl 



+ 



gjy 

8tt 2 



/2 



( 7 2 + 2 V /^ 2 -^-<7 2 )+^ + C. 



(50) 



Next, we have to find the minimum point a of this effective potential, i.e. to solve the stationarity 
equation dV^jda 2 = 0, which gives for a the equation 



3A 2 3a, 



3// 

"27 2 2tt 2 4tt 2 



4tt 2 8 2A 2 16tt 2 ^ 



2/i' 



^^' 2 -gH 4G 



+ 77^ + O(^ 2 )=0. 
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It is seen that, contrary to the case with the pion condensate, there are no large terms like 
gH\og(gH/a 2 ) in the above equation, and hence, for large terms with A s , it has a nontrivial solution 
cr ^ only if g > 1, which is the same condition as for the zero field case (see, e.g., [HI), 

a = VdAexp - h) > (51) 

where C is a certain numerical constant of order unity. 

Comparing the results of the particular considerations a) and b), the following preliminary con- 
clusions can be made. It is clear that at very small values of the coupling constant G, i.e. at g < 1, 
and for // > 0, both the isotopic and the chiral symmetry of the model (0) are not broken down 
at H = 0, and quarks are massless particles (recall, we consider only the case with zero current 
quark mass). However, if an external chromomagnetic field H is present (and it is rather strong in 
our analytical consideration), then the effective potential (J23J) of the system acquires a nontrivial 
global minimum point which lies on the 7Ti-axis and has the form (a = 0, 7Ti = ttq ^ 0), where 
7To is the solution of the stationarity equation (|43|). This means that in this case a nonzero pion 
condensate (q^^Tiq) = —tiq/2G is generated by an external chromomagnetic field (chromomagnetic 
catalysis of the pion condensation), and the isotopic symmetry C// 3 (l) of the model is spontaneously 
broken down (see the text after Moreover, since the expression for the effective potential in this 
case resembles the effective potential in the 2-dimensional Gross - Neveu model, one can say that 
this effect, similar to the magnetic catalysis phenomenon, is provided by the dimensional reduction 
mechanism. 

Let us present some numerical calculations in the next section, in order to support this conclusion. 



III. NUMERICAL RESULTS 



Consider first the case of a zero external field F^ v = 0, and flavor symmetric medium // = 0. 
Moreover, the temperature and bare quark mass m are taken to be zero throughout the present 
section. The numerical analysis of the behavior of the effective potential (|25|). which is just the 
thermodynamic potential Q for T = 0, was performed for the value of the coupling constant G = 5.01 
GeV -2 , taken from j24[ 25[, and for the quark chemical potential /i = (see Fig. Q). In this case the 



quark energy spectrum is symmetric with respect to a and ir%, i.e. the function V e ff(a, tti) depends 
only on the single variable \/ a 2 + irf. Hence we draw the effective potential as a function of only 
one variable, e.g., a. In this figure, the picture on the left corresponds to the value of the cutoff 
parameter A = 0.165 GeV (in this case g = 0.08), whereas the right picture corresponds to A = 0.65 
GeV (g = 1.29). The figure demonstrates the general feature of the effective potential at y! = and 
gH = 0. Namely, if the effective coupling constant g is sufficiently small, g < 1, then the global 
minimum point of the effective potential is at the origin, and the chiral symmetry is not broken (left 
picture). However, if g > 1, the effective potential has a nontrivial global minimum, and the chiral 
symmetry is spontaneously broken down (right picture). The dependence of V e g(a, tti) on the two 
variables a and it\ (g = 1.29) is depicted in Fig. |2]for fi = 0. The picture is evidently symmetric, as 
it should be in this flavor symmetric case. 

It is also necessary to note that if g > 1, then at g H = and |i'>0a nonzero pion condensate 
appears in the model (JTJ) [i^, 25l ■ (In this case, the isotopic symmetry £/j 3 (l) is broken down). 
However, if g < 1, then at gH = and arbitrary values of including zero, the symmetry of the 
NJL model (0) remains intact. 

Now, in order to present numerical arguments in favour of the chromomagnetic catalysis of pion 
condensation, we will study the effective potential V e s(a, 7Ti) (j2*B]) at //' > and suppose that 
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gH = 0.5 GeV 2 (this value of gH mimicks the nonzero QCD gluon condensate [13 )■ Since 
the current quark mass m is equal to zero, the nontrivial stationary points of the function (|25p 
lie either on the a or on 7Ti axis (see the remark following formula (|15|)). Thus, in order to find 
the global minimum point of the effective potential (J2*5j) . it is enough to consider the behavior of 
V e g(a, 7Ti) along the coordinate a and 7ri-axis only. The corresponding curves are presented in Fig. El 
for fi' = 0.15 GeV, and physical values of G = 5.01 GeV -2 and A = 0.65 GeV (in this case 
g = 1.29). They demonstrate that, although a nontrivial stationary point of the thermodynamic 
potential Q = V e g at T = does exist on the cr-axis in a chromomagnetic field gH (and this is 
just a saddle point in this case), the value of the effective potential at the other stationary point, 
a = 0, TV i = 7i 7^ 0, is slightly deeper. Our calculations indicate that with growing // the differ- 
ence between the minima of the curves for V e e(iri) and V e s(a) increases, and the very minima get 
deeper. The comparison of the results obtained with the known results at zero background field, 
demonstrated that the background field deepens the minimum for the pion condensate. Moreover, 
we formally considered the interesting case of a smaller value of A, corresponding to g < 1. The 
calculated effective potential then has a minimum at a nontrivial value of 7r for weak coupling of 
quarks for g < 1, as well. The corresponding curves are presented in Figs. EJ |5]for different values 
of fjf = 0.1 GeV, // = 0.15 GeV, respectively, and G = 5.01 GeV" 2 and A = 0.46 GeV (in this case 
g = 0.64). 4 Thus, only at this point, a = 0, 7Ti = ttq ^ 0, the global minimum of the thermodynamic 
potential (|2*3jl is observed, and the nonzero pion condensate (q^Tiq) = —■kq/2G is generated by the 
external chromomagnetic field in the model. 




sigma 



FIG. 1: Thermodynamic potential Q, = V e ff(a, 0) at T — 0, // = 0, /i = and G = 5.01 GeV 2 as a function of a (in GeV-units) 
for A = 0.165 GeV - left picture, and A = 0.65 GeV - right picture. 



IV. CONCLUSIONS 

We considered, in the framework of an NJL model, the effects of quark and pion condensation 
in dense quark matter with and without flavor asymmetry under the influence of an external chro- 
momagnetic field modelling the gluon condensate (F^ u F atJ-u ). The general conclusion is that the 
presence of the chromomagnetic field catalyses the effect of quark or pion condensation, depend- 
ing on the values of the chemical potential jj! . Indeed, when there is no chromomagnetic field, no 



4 In contrast with the numerical results, in our analytical investigations (see the previous section), the stationary point of the effective 
potential was not observed on the tr-axis at g < 1. This fact can be explained due to a different domain for the variable fi' used in these 
cases. Indeed, it was supposed in section II. B that fi' = O(A), whereas for the numerical consideration, we used fi' <g A. 
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FIG. 2: Thermodynamic potential f2 = V e g(a, 7Ti) at T = 0, = 0, // = as a function of a, 7Ti for A p = 0.65, /i = 0. 
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FIG. 3: Thermodynamic potential Q = Veff(ff,iri) at T = = 0, t/i? = 0.5 GeV 2 , A = 0.65 GeV, g = 1.29, fj,' = 0.15 GeV as 
a function of a (in GeV-units) at 7Ti = - left picture, and of tti (in GeV) at <r = - right picture. 




FIG. 4: Thermodynamic potential fi = Kff(cr,7ri) at T = 0, /x = 0, gH = 0.5 GeV 2 , A = 0.46 GeV, g = 0.64, fx' = 0.1 GeV as 
a function of a (in GeV-units) at 7Ti = - left picture, and of 7Ti (in GeV) at a = - right picture. 
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FIG. 5: Thermodynamic potential Q = V eB (a, tti) at T = 0, /x = 0, gH = 0.5 GeV 2 , A = 0.46 GeV, g = 0.64, fx' = 0.15 GeV as 
a function of a (in GeV-units) at 7Ti = - left picture, and of iri (in GeV) at a = - right picture. 



condensation of quarks or pions takes place at low values of the effective coupling (J3l)j) of quarks 
g < 1. However, when the chromomagnetic background field is present, the pion (at fi' > 0) or 
quark condesation (see the case // = of Section II) take place even for small values of g < 1 (recall 
that in our consideration the current quark mass m was taken to be zero). If g > 1, the role of 
the external chromomagnetic field H is to enhance the effect of pion condensation that appeares in 
the model at // > already at H = 0. For further investigations, it would be interesting to study 
the influence of an external chromomagnetic field (the gluon condensate) on the pion condensation 
phenomenon at nonzero current quark mass m as well as to consider the chromomagnetic catalysis 
effect with another combinations of a background field, different from (|20|). 
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